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Abstract 

We study the electronic transport in a graphene-based superconductor-normal(graphene)- 
superconductor (SNS) junction by use of the Dirac-Bogohubov-de Gennes equation. We consider 
the properties of tunnehng conductance through an undoped strip of graphene with heavily doped 
superconducting electrodes in the dirty limit l^ef <^ L <^ ^. We find that spectrum of Andreev 
bound states are modified in the presence of single localized defect in the bulk. The minimum 
tunneling conductance remains the same and this result doesn't depend on the actual location of 
the imperfection. 
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I. INTRODUCTION 



Graphene, namely a monolayer of graphite, is formed by carbon atoms on a two- 
dimensional honeycomb lattice. In graphene, due to its unique band structure whose valence 
and conductance bands touch at two inequivalent Dirac points (often referred to as K and 
K') of the Brillouin zone, the electrons around the Fermi level obey the massless relativis- 
tic Dirac equation which results in the linear energy dispersion relation. Recent exciting 
developments in transport experiments on graphene has stimulated theoretical studies of 
superconductivity phenomena in this material, which has been recently fabricated HQ. A 
number of unusual features of superconducting state have been predicted, which are closely 



related to the Dirac-like spectrum of normal state excitation [3|, 14] • In particular, the uncon- 
ventional normal electron dispertion has been shown to result in a nontrivial modification of 
Andreev reflection and Andreev bound states in Josephson junctions with supercomducting 



graphene electrodes 



Other interesting consequences of the existence of Dirac-like quasiparticles can be under- 
stood by studying superconductivity in graphene js- 12] . It has been suggested that super- 
conductivity can be induced in graphene layer in the presence of a superconducting electrode 



near it via proximity effect 13l-ll5|. 



In this work, we study Josephson effect and find bound state in graphene for tunnel- 
ing SNS junction with the presence of a single localized defect 16|. In this study, we shall 
concentrate on SNS junction with normal region thickness L ^ ^, where ^ is the super- 
conducting coherence length, and width W which has an applied gate voltage U across the 



normal region 



13, 



18j. In the frame of the dirty limit Idef L C, cosidered by Kulik 



and Omelyanchuk we investigate tunneling condactance in SNS junction with presense a 
single localized defect and find that Andreev levels are modified, the minimum tunneling 



conductance remains the same 



19|421|. 
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FIG. 1: A graphene undoped ribbon is contacted by two superconducting leads. The charge carriers 
tunnel from one lead to another via multiple tunneling states formed in the graphene strip. An 
defect placed inside the strip. 

II. TUNNELING CONDUCTANCE OF THE GRAPHENE SUPERCONDUC- 
TOR/NORMAL/SUPERCONDUCTOR JUNCTION WITH A SINGLE LOCAL- 
IZED DEFECT 

We consider a SNS junction with a single localized defect which is involved in a graphene 
sheet of width W lying in the x — y plane extends from x = — L/2 to x = L/2 while the 
superconducting region occupies |x| > L/2 (see Fig. 1). The SNS junction can then be 



described by the Dirac-Bogoliubov-de-Gennes (DBdG) equations[22 

'4's = e^, 



Hs-Ef + U a 

A* Ep-U-H, 



Here , ips = {i'As,4'Bs,4'Asy ~'^Bs)y V' = (""i; '?^2, "^^i, "^^2) are the 4 component wave functions 
for the electron and hole spinors, the index s denote K or K' for electrons or holes near K 
and K' points, 1 takes values K{K') for s = K(K'), Ep denotes the Fermi energi, A and B 
denote the two inequivalent sites in the hexagonal lattice of graphene, and the Hamiltonian 
Hg is given by 

Hs = -ihvF[(Jxdx + sgn{s)aydy\ (1) 

In Eq. [H vf denotes the Fermi velocity of the quasiparticles in graphene and sgn{s) takes 
values ± for s = K(K'). The 2x2 Pauli matrices cij act on the sublattice index. The 
excitation energy e > is measured relative to the Fermi level(set at zero). The electrostatic 
potential U and pair patential A have step function profiles, as in the case of a semiconductor 
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two-dimensional electron gas 23N25| . 



U{x) 



-U, if X < -L/2, 
0, if |a;| < L/2, 
-U, if X > L/2. 



Aoexp(z0/2), ifx<-L/2, 
A(x) = < 0, if < L/2, 

Aoexp(-i0/2), ifx>L/2. 
The reduction of the order parameter A(a;) in the superconducting region on approaching 
the SN interface is neglected; i.e., we approximate parameter A(x) as we've done it above. 



As discussed by Likharev [26[], this approximation is justified if the weak link has length 
and width much smaller than C,- There is no lattice mismatch at the NS interface, so the 
honeycomb lattice of graphene is unperturbed at the boundary, the interface is smooth and 
impurity free. 

Solving the DBdG equations, we gain the wave-functions in the superconducting and 
the normal regions. In region I {1 1 1), for the DBdG quasiparticles moving along the ±x 
direction with a transverse momentum ky = q and energy e, the wave-functions are given by 

/ exp(— zm/3) \ 
exp(z7 — imP) 



\Ef = exp{iqy + ikgX + nmx) 



exp(— zm0/2) 
y exp(i7 — im(f)/2) J 



expiiqy — ikgX + nmx) 



( exp(im/3) \ 
exp(— Z7 + im0) 
exp(— im0/2) 
y exp(— Z7 — im(f)/2) J 
The parameters (3, 7, k^, k are defined by /3 = arccos(e/Ao), 7 = arcsin[^t>i?g/(f/o + Ep)], 
ks = A/(f/o + Epf/iflVFY - q^, K = (f/o + EF)Aosm{(3)/{hhlks) and m = ± denotes 
region I {1 1 1), m = + for / and m = — for /// correspondent ly. Further we assumed that 
the Fermi wave length in the superconducting region much smaller than the wave length 
Xp in the normal region and Uq ^ Ep,e. Since |g| < Ep/hvp, this regime of a heavely 
doped superconductor corresponds to the limits 7 — )■ 0, /c^ — t- Uo/hvp, n — )• {Aq/Hvp) sin(/3). 



FIG. 2: The normal region II is devided by three areas a,b,c. The DBdG equations are solved 
for area b and determined by the "infinite mass" boundary conditions induced by V{r) — > +00 in 
a and c areas correspondently. 



The region // consists of three areas: a,b,c (see Fig. 2). We solve the DBdG equations 
for the area 6, while area c is the area where placed defect and area a would be extended 
and matched with superconducting regions. The two valleys s± decouple, and we can solve 
equations separately for each valley, Hgi's = (e + sEp)ips, Hs = Hq + sV{r)az. The term 
proportional to cr^ in Hamiltonian is a mass term confining the Dirac electrons in the area 
b. Rewrite the Hamiltonian in the cylindrical coordinates and since Hg commutes with 



Jz = h + ^<^z, its electron-eigenspinors ipe are eigenstates of [27l |. 

/ exp{id{n-l/2)a)Jd(n-i/2){k{e)r) 

y exp{id{n+ l/2)a)Jd(n+i/2)iH<^)r) 
with eigenvalues n, where n is a half-odd integer, n = d^,d^,... and Jd(n-i/2){k{e)r) is 
the Bessel function of {n — 1/2) order. In the x — y plane d denotes moving direction of 
correspondent quasiparticle, ci = + for the quasiparticle moving toward x = L/2 and d = — 
for the quasiparticle moving toward x = —L/2 direction correspondentely. Further we are 



interested in to find zero energy states[15|]- In this case the DBdG equations posses a general 



symmetry with respect to the change in the sign of energy 

e — )■ — e, idyU* — > v, icfyV* — )■ — m, (2) 

where we denote u = (ui,M2) and v = (^1,^2). Thus, for a set of zero modes {ui,Vi) 
enumerated by a certain index i we should have, 

Vi = iayU*, Ui = —idyV*. (3) 



FIG. 3: Scheme of graphene SNS junction with a single impurity which is placed inside the strip 
of the region II. 



In the same manner as for electron hole-spinors have form 

-exp(zci(n + l/2)a')Jd(n+i/2)ik'{e)r) 
ex.p{id{n - l/2)a')Jd(n-i/2)ik'{e)r) 



with the definitions 

a{e) = a.Tcsm[hvFq/{e + Ep)], a (e) = a.Tcsm[hvFq/{e — Ep)], (4) 

k{e) = {hvF)~^{e + Ep) cos(a), k (e) = (/ifp')~^(e — Ep) cos(a;), (5) 

The angle a G (— 7r/2,7r/2) is the angle of incidence of the electron (having longitudinal 
wave vector fc), and a is the reflection angle of the hole (having longitudinal wave vector k') 
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29| . To obtain an analytical approximation of the spectrum, we use the asymptotic form 
of the Bessel functions for large r. This indeed is the desired limit as rk{e) ~ Tdefk{e) oc 
TdeflL <^ 1, where r^e/ is the defect radius (the radius of the area c) and determine for 
all eigenvalues n = d^. In this limit we impose the "infinite mass" boundary conditions at 
y = 0, W, for which Qn = {n + 1/2)71 /W in the area h with V{r) — > +oo in the a and c areas 
consequently . Half-odd integer values n reflect the vr Berry's phase of closed size of a single 
localized defect in graphene. 

To obtain the subgap (e < Aq) Andreev bound states, we now impose the boundary 
conditions at the graphene. The wave-functions in the superconducting and normal regions 



can be constructed as 



= aiipf + hi) J , -^^iii = a2ipt + b2tpi , (6) 
= a^ii + briT + cV/Y + . (7) 
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where ai (61), 02(62) are the amphtudes of right and left moving DBdG quasiparticles in 
region I {1 1 1) and a{b) and c{d) are the amphtudes of right (left) moving electrons and 



holes, respectively, in the normal region js]. These wave functions must satisfy the boundary 
conditions, 

^l\x=-L/2 = ^//U=-L/2, ^Il\x=L/2 = ^IIl\x=L/2- (8) 

Since the wave vector ky parallel to the NS interface and different wave vectors in the y- 

direction are not coupled, we may solve the problem for a given ky = q and we can consider 

each transverse mode separately. To leading order in the small parameter AqL/Hvf we may 

substitute a{a ) ^ a(0), k{e){k (e)) ^ A;(0). After some algebra we obtain equation 

, ^, f sin^(A;L) — cos^(/cL) sin^(a) 1 , ^, f cos(A;L) sin(A;L) sin(a) 1 , ,, , , 

-°-P^> I J(kL) cos» - 1 I - -'■'P'^' I J^-(kL) Lh.) - 1 I - 

Eq. [9] differs from equation obtained for SNS junction without a single defect. Eliminating 

second term in Eq. [9] we could immediately yield reduction of the equation and it earns a 



essential form for weak SNS junctions 
mode, 



15|. The solution of Eq. [9] is a single bound state per 



= Ao^^^-^i-2CA^ -B^ + B^l-ACA\C + l)) (10) 

where 

^ sin^(/c„L) — cos^(A;„L) sin^(Q;) ^ cos{knL) sin(/c„L) sin(a) 
cos^(/c„L) cos^(q;) — 1 ' cos^(/c„L) cos^(q;) — 1 

C = i + 1(1 - .in^(|)). r„ ^ ,,7;f -^1;^' °' (12) 

2 r„ 2 2 cos^(A;„L) cos^(a) — cos(2A;„L) 

We don't have a simple analytic expression for the 0-dependance but we obtained modified 
Andreev levels with the presense a single localized defect in the bulk. The conductance of the 
graphene strip is expressed through the transmission probability by the Landauer formula, 

G = go^Tn, gQ = Ae^/h, (13) 

n=0 

where n(/x) ^ 1 is given by n{^)=hii{knW / + 1/2). Substitution transmission probability 
into Eq. [13] gives the conductance versus Fermi energy (see in Fig. 3, Fig. 4). The result 



for the minimal conductivity agrees with other calculations 30l-l32l| , which start from an un- 
bounded disordered system and then take the limit of infinite mean free path /. There is no 
geometry dependence if the limits are taken in that order. 




FIG. 4: Tunneling conductance of the graphene SNS junctions with a single localized defect versus 
Fermi energy, calculated from Eq. [T2l The tunneling conductance exhibits oscillatory behavior. 

III. SUMMARY 

In the conclusion, we have shown that tunneling conductance of the graphene SNS junc- 
tion with a single localized defect has a nonzero minimal value if the Fermi level is tuned 
to the point of zero carrier concentration. We have demonstrated that the tunneling conduc- 
tance exhibits oscillatory behavior. Andreev levels are modified, the minimum tunneling 
conductance remains the same and this result doesn't depend on the actual location of the 
imperfection. The consideration of tunneling conductance of SNS junctions with multiple 
defects is the subject of forthcoming article. 
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